  function res=mat4bc(ya,yb,lambda)
  res=[ya(2)
        yb(2)
        ya(1)-1];
